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SUMMARY 


A combined analytical and experimental determination is made of the 
coupled natural frequencies and mode shapes in the longitudinal plane of 
symmetry for a dynamic model of a single-rotor helicopter. The analyti- 
cal phase is worked out on the basis of a seven -degree -of -freedom system 
combining elastic deflections of the rotor blades, rotor shaft, pylon, 
and fuselage. The calculated coupled frequencies are first compared with 
calculated uncoupled frequencies to show the general effects of coupling 
and then with measured coupled frequencies to determine the extent to 
which the coupled frequencies can be calculated. The coupled mode shapes 
are also calculated and were observed visually with stroboscopic lights 
during the tests. 

A comparison of the coupled and uncoupled natural frequencies shows 
that significant differences exist between these frequencies for some of 
the modes. Good agreement is obtained between the measured and calcu- 
lated values for the coupled natural frequencies and mode shapes . The 
results show that the coupled natural frequencies and mode shapes can be 
determined by the analytical procedure presented herein with sufficient 
accuracy if the mass and stiffness distributions of the various compo- 
nents of the helicopter are known. 


INTRODUCTION 


It has become a well substantiated truth (for example, refs. 1 to 6) 
that helicopter rotors, when moving in forward flight, develop periodic 
aerodynamic loads which possess harmonic components having frequencies 
equal to all lower integral multiples of the rotational frequency of the 
rotor (the first 10 usually being of fundamental significance). It is 
also well known that the relative magnitudes of the harmonic components 
of the aerodynamic loading vary with rotor configuration and are strongly 
dependent on the helicopter flight velocity. These periodic aerodynamic 
forces produce undesirable vibrations in the helicopter which, in 
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conformity with the theory of vibrations, are dependent on the magnitude 
of the exciting force, the inherent or applied damping of the structure, 
and the nearness to resonance of the frequencies of the exciting forces 
with the natural frequencies. 

The aerodynamic input forces and structural damping forces are con- 
trollable only to a limited extent and, thus, for practical purposes, 
the .best way to avoid undesirable vibrations is to avoid conditions of 
resonance. In order to do this in a systematic and efficient manner, it 
is necessary to be able to predict the natural frequencies of the heli- 
copter in flight. The frequencies referred to here are the coupled 
frequencies of the entire helicopter and not just the frequencies of the 
individual components . 

A combined flight test and analytical investigation for a tandem 
helicopter is given in reference 7> wherein the natural frequencies, 
mode shapes, and vibration amplitudes in the longitudinal plane of sym- 
metry are presented. The vibration amplitudes are presented for both in- 
flight aerodynamic forces and applied eccentric mass shaker forces. In 
the analytical investigation of the coupled natural frequencies therein, 
an attempt was made to select the more significant variables appropriate 
to that configuration in order to study the overall effects of coupling. 
Correlation between experimental and analytically predicted frequencies 
was good, but because the helicopter was not disassembled to obtain 
accurate evaluations of the mass and stiffness distributions of its com- 
ponents and inasmuch as the coupled mode shapes were not measured during 
flight tests, no precise evaluation of the adequacy of the assumptions in 
the theory could be made. Therefore, the decision was made to extend the 
investigation by means of dynamic-model studies wherein the structural 
properties of the components and the experimental coupled natural fre- 
quencies and mode shapes could be measured with sufficient accuracy to 
permit an evaluation of the theoretical approach. 

This report presents the results of a combined experimental and 
analytical study of the coupled natural frequencies and mode shapes of a 
dynamic model. The specific purpose is to demonstrate that the general 
effects of coupling as outlined in reference 7 do exist, and that it is 
possible to predict, with reasonable accuracy, the coupled natural fre- 
quencies and mode shapes of the helicopter if the structural properties 
of the components are known. Although the experimental portion of the 
study was made under hovering conditions, the results are appropriate 
for any flight condition because the change in coupled frequencies and 
mode shapes with changes in flight conditions are believed to be 
negligible. In the analytical portion of this study, the theory of 
reference 7 is extended to include lateral-flexibility characteristics 
of the pylon and shaft. 
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In the organization of this report, the description of the model is 
first given in some detail. The reader who is interested primarily in 
the analysis may refer directly to the section entitled "Analytical 
Determination of Modes and Frequencies." 


SYMBOLS 




a~L y 


c i 

d o' d o> d l 

e 

(El) f 

( EI )p 

(El) r 

F 

f 

K s 

K x 


k l ,k 2 


Mp 


coefficients of frequency determinant 
coefficients of rotor-blade modal deflections 

coefficient of rotor-shaft modal deflections 
coefficient of pylon modal deflections 
coefficient of fuselage modal deflections 

flapping hinge offset, in. 

structural stiffness of fuselage, lb-in. ^ 

p 

structural stiffness of pylon, lb-in. 

p 

structural stiffness of rotor blades, lb-in. 
length of fuselage, in. 

longitudinal position of any fuselage element, in. 

spring constant for shaft motions relative to the pylon, 
lb/in. 

spring constant of mounts for vertical translation of 
fuselage, lb/in. 

spring constant of mounts for pitching motion of 
fuselage, lb-in. /radian 

Southwell coefficients for flapwise blade bending 
mass of fuselage, lb-sec^/in. 
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M r 


mass of pylon, lb-sec^/in. 


Mq 

“R 

M d 


M j = M 
sd s 


M M 
1’ 2’ ’ 

TTL^* 

m 

P 


m v 


n l ,n 2 . 


mass of rotor shaft, lb-sec^/in. 

mass of rotor blades, lb-sec^/in. 

mass of dampers in upper pylon frame, lb-sec 2 /in. 

mass of springs in upper pylon frame, lb-sec^/in. 

+ M d 

. . apparent mass terms in frequency determinant 

mass per unit length of fuselage, Ib-sec^/in.^ 
mass per unit length of pylon, lb-sec^/in.^ 
mass per unit length of rotor shaft, lb-sec^/in.^ 
mass per unit length of rotor blades, lb-sec^/in.^ 
. . coefficients in frequency determinant 


n 

P 

P 

Q 

q 

R 

r 

T 

T 


number of blades 

vertical position of upper pylon frame, in. 
vertical position of any pylon element, in. 
vertical position of rotor, in. 

vertical position of any rotor shaft element, in. 
rotor radius, in. 

radial position of any blade element, in. 

centrifugal force on the blade at any radial station, lb 

centrifugal force on the blade at any radial station 
divided by square of rotor speed, lb-sec^ 
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u 


u l> u 2 

V 

v o 

W 1 

x o>*o' x l 


kinetic energy, lb-in. 

uncoupled mode shapes for rotor blades 

potential energy, lb-in. 
uncoupled mode shape for a rotor shaft 
uncoupled mode shape for pylon 
uncoupled mode shapes for fuselage 


n 

O) 

“( 1 )'*( 2 )' * 
o>i 

cx>2 


effective mass ratio 

rotor speed, radians/sec 

natural frequency of coupled modes, radians/sec 

natural frequencies of first coupled mode, second 
coupled mode, . . . 

first natural uncoupled nonrotating flapwise bending 
frequency of rotor blades, radians/sec 

second natural uncoupled nonrotating flapwise bending 
frequency of rotor blades, radians/sec 

first natural uncoupled vertical bending frequency of 
fuselage , radians /sec 

first natural uncoupled bending frequency of pylon, 
radians/ sec 


Subscripts : 

a intersection of rotor shaft and rotor-disk plane 

p intersection of rotor shaft and upper pylon frame 

y intersection of rotor shaft and horizontal plane of 

center of fuselage 

Dots denote differentiation with respect to time, and primes denote 
differentiation with respect to length or space. 
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APPARATUS AND TESTS 


General Description of Model 

A general view of the model used in this investigation is shown in 
figure 1 , details of the various components are shown in figures 2 to 6 , 
and a sketch of the test configuration is shown in figure 7* The model 
was designed to possess the approximate dynamic properties of a l/ 8 -scale 
model of an existing helicopter. In this respect, the flapwise bending 
stiffness and mass distribution of the blades and the bending stiffness 
and mass distribution of the fuselage and pylon were scaled to yield a 
frequency spectrum representative of existing helicopters. 

The model was designed for general dynamics studies, including the 
one discussed in the present report, and contains many provisions for 
variation of the significant parameters such as mass and stiffness 
distributions . 


Description of Components 

Rotor . - The rotor is a four-blade configuration with fully articu- 
lated blades. (See figs. 1 , 2, and 4.) Motion of the blades about the 
drag hinges is damped by small viscous shear-type dampers (see figs. 3 , 

4, and 5) to avoid ground-resonance instabilities. Positive stops are 
also incorporated to limit blade motion in the plane of rotation to ± 12 °. 
The collective pitch angles of the blades are adjusted individually by 
rotation of the blades in pitch at the blade root clamp. The pitch axis 
is located at 25 percent chord. 

The rotor blades are of rectangular plan form and have a chord of 
2.062 inches. The rotor has a diameter of 66 inches and a solidity of 
0.08. The airfoil section is an NACA 0015 section. The blade spanwise 
weight distribution outboard of the root attachment is uniform with a 
chordwise center-of -gravity location at 32 percent chord. The blade 
structure consists of an aluminum leading-edge spar (see fig. 6) with 
balsa laminated to the spar to form the airfoil section. The flapwise 
bending stiffness of the rotor blades was scaled to be representative 
of existing helicopter rotor blades, whereas the torsional rigidity was 
maintained rather high to minimize the possibility of coupling between 
the flapwise bending and blade torsion. A trim tab attached to the 
trailing edge of the blade extending from station 30*5 to 31.5 and con- 
sisting of a small-gage wire frame covered with linen allows adjustment 
for blade tracking. Provisions for adjustment of the angular spacing 
between the blades and mass balance of the rotor system are obtained by 
the addition of weights to a balance rod shown in figure 4 . 
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Pylon and rotor shaft .- The details of the pylon and rotor shaft 
are shown in figures 2 and 3- The pylon consists of a lower frame and 
upper frame which are connected by four vertical rods. The lower frame 
is rigidly attached to the fuselage . The upper pylon frame and the rotor 
shaft bearing housing are connected by four springs and dampers as shown 
in figures 2 and 3 * The rotor shaft is connected at the top of the gear 
box (in the plane of the lower pylon frame) by means of a universal joint 
which permits the upper portion of the shaft to translate relative to the 
pylon. Thus , the flexibility of the shaft relative to the pylon is con- 
trolled by the four springs in the upper pylon frame. The flexibility of 
the entire pylon-shaft system relative to the fuselage is controlled by 
the stiffness of the four vertical pylon rods. This combination of 
springs and vertical rods provides a convenient method for variation oi 
the natural frequencies of the pylon-shaft assembly. 

The gear box at the lower end of the pylon is rigidly attached to 
the fuselage and transmits the torque from the horizontal shaft in the 
fuselage to the rotor shaft and, thence, to the rotor. 

Fuselage and fuselage suspension .- The fuselage, shown in figure 1, 
was of box cross section and was constructed of two welded magnesium angles 
to provide the desired stiffness. The desired mass distribution is then 
obtained by the addition of mass as required. In simulating the fuselage, 
an attempt was made to obtain desired dynamic characteristics only, the 
geometry of the fuselage being secondary inasmuch as the rotor thrust 
was small and the aerodynamic effects on the natural frequencies and 
mode shapes negligible. 

The fuselage was mounted on a very stiff cantilever beam, as shown 
in the sketch of figure 7 , by means of two supports (front support and 
rear support) which permitted the model freedom of vertical translation 
and longitudinal pitching but provided a high relative restraint to other 
rigid-body motions. These supports were located at the approximate posi- 
tions of the fuselage node points for fuselage first elastic free -free 
bending. The objective In using this type of support system was to 
duplicate as closely as possible the characteristics of a free -free 
system in the longitudinal plane of symmetry. During the model tests 
the configuration of the front support was found to be more desirable 
than the rear support (as shown in fig. l) , and subsequently the rear 
support was changed to be identical to the front support. The data pre- 
sented in this paper were obtained with the latter support configuration 
(as shown in fig. 7)* 

Dampers . - The small dampers used to dampen motion of the blades 
about the drag hinge and the motion of the rotor shaft at the upper 
suspension plate of the pylon are of the shear type. The damper consists 
of a rod which slides through neoprene tubing. (See fig. 5*) A 
detailed description of the characteristics of this type, of damper is given 
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in reference 8. Silicone oil is used as the damping fluid and the 
dampers provide a damping force which varies approximately as the square 
root of the velocity. The blade dampers also contain stops to limit 
rotor-blade -lag angles. 

Model power supply .- The rotor system is driven by a compressed-air 
motor mounted in the center of the fuselage. A static pressure of 
150 lb/sq in. is regulated by a gate control valve from which the air 
passes through a short flexible line to isolate vibration sources to the 
airline. The power output of the air motor is transmitted through a 
horizontal drive shaft to a commercial right-angle gear box which drives 
the rotor shaft. 

Instrumentation . - Instrumentation of the model was accomplished by 
the use of accelerometers and electrical resistance -wire strain gages. 
The rotor-blade strains were indicated by flapwise bending and torsion 

strain gages located at £ = O. 38 . Strain gages mounted on small canti- 
lever beams were used to measure pylon motions relative to the fuselage 
and rotor-shaft motions relative to the pylon. (See fig. 3*) Vertical 
translation and pitching motions of the fuselage were measured by accel- 
erometers mounted on the fuselage and strain gages located at fuselage 

stations of - = 0.33 and 0.66 indicated fuselage vertical bending. 

F 

The strain-gage signals from the rotor are transmitted to the 
recording equipment through a slipring-brush assembly mounted on the 
rotor shaft as shown in figure 2. 

An instantaneous value of the rotor speed was obtained by means of 
a commercial rotational -speed-indicating device which received signals 
from the small multipole generator operated from the gear box and shown 
just below the pylon in figure 1. A permanent record of the rotor speed 
was obtained from the signal of a one -per-re volution indicator mounted 
on the generator drive shaft. 

Stroboscopic lighting was employed throughout the tests to obtain 
qualitative definition of the mode shapes associated with the various 
coupled natural frequencies . 

All data were recorded on oscillograph records. 


Evaluation of Component Structural Properties 

Whenever possible in the method of analysis presented in this report, 
the equations are set up in a manner that permits the inclusion of the struc- 
tural properties in terms of the spring constants, natural frequencies, 
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and mass distributions of the uncoupled components (blades, shaft, pylon, 
fuselage, and so forth). These component properties can, of course, be 
determined either experimentally or analytically, and in some instances, 
both procedures were used. However, greater emphasis was given to the 
exDerimental approach because it is felt that a similar procedure would 
lead to a more satisfactory evaluation of a full-scale helicopter. For 
example, the potential energy associated with the deflection of a rotor 
blade in a given mode shape is more readily determined by measuring the 
natural frequency and calculating the generalized mass for the particular 
mode than by integrating the product of the structural stiffness and the 
square of the blade curvature. In some instances, however, the experi- 
mental values were checked by the analytical methods presented in refer- 
ences 9, 10, and 11 to assure that these methods could be used to obtain 
sufficiently accurate results for preliminary calculations before com- 
ponents are available for testing or in the event that experimental test 
equipment is not available for component testing. 

Prior to the assembly of the model, the natural frequencies of the 
individual components were determined experimentally . This was accom- 
plished by mounting the component in such a manner as to simulate the 
assembled configuration (fixed-free condition for pylon, pin-end condi- 
tion for blades, and so forth). The natural frequencies were then 
obtained by exciting the components with suitable shakers and noting the 
frequencies, as obtained from oscillograph records, where the deflection 
became a maximum. 


Measurement of Coupled Natural Frequencies 
and Mode Shapes 

After assembly of the model, the I-beam shown in figure 1 was 
mounted in a horizontal position on a rigid test bed in cantilever 
fashion as shown in figure 7. At a selected rotor speed, the frequency 
of the shaker output force was slowly varied from near zero up to the 
desired maximum value of approximately 300 radians per second. As the 
frequency was varied, the response of the various components was 
measured and recorded on oscillograph records. A simultaneous visual 
observation of the model was obtained by means of stroboscopic lights 
to obtain a qualitative value of the coupled natural frequencies as 
resonance conditions arose and to observe the corresponding coupled mode 
shapes. The rotor speed was then changed and the process repeated in 
order to obtain the variation of the coupled natural frequencies with 
rotor speed. 

The coupled natural frequencies for each rotor speed were obtained 
from visual observations of the oscillograph records. A condition of 
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resonance vas considered to exist at the frequencies where the response 
of the system became a maximum. These frequencies were then checked, 
wherever possible, against the values which produced maximum model 
deflections as determined by visual observations of the model under 
stroboscopic lighting. These visual observations were quite good for 
those modes involving large relative amplitudes of fuselage, pylon, and 
shaft motions but were of little use for modes involving large relative 
blade motions . 


ANALYTICAL DETERMINATION OF MODES AND FREQUENCIES 


General Considerations 

The analysis of the coupled frequencies and mode shapes of the 
model employs the Lagrange energy -equation approach of reference 7. The 
analysis consists of writing the kinetic and potential energies of the 
structure in terms of the displacements of individual components . 

Although dampers were installed on the model to avoid ground -resonance 
instabilities, the effect of damping on the coupled frequencies is 
assumed to be of secondary importance. The calculation of the natural 
coupled frequencies is limited to the following degrees of freedom in 
the longitudinal plane of symmetry in order to reduce the complexity of 
the derivation: vertical translation, longitudinal pitching and struc- 

tural bending of the fuselage, first and second flapwise bending of 
rotor blades, structural cantilever bending of the pylon, and rotation 
of the rotor shaft about the universal joint. In addition, an allowance 
is made for blade chordwise motions by assuming that the effective mass 
of the rotor blades in fore-and-aft translations of the rotor hub is qM^ 

where q is generally less than unity. In the calculations of this 
paper, the effective mass ratio q is assumed to be 0.75- 

Substitution of the energy equations into Lagrange* s dynamic 
equations for free vibration yields the equations of motion from which 
the coupled frequencies and mode shapes are determined. 


Energy Equations 

The kinetic and potential energies are written in terms of the 
deflections of the components and their derivatives as follows: 
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Kinetic energy U 


U = a 
2 


'R 


m T 


UO 


(il + X 7 ) 2 + T](v a ) 2 


ro. 

dr + ^ m 

^ Jo 


^v 2 + dq 


1 

2 


f F 

mp^v 2 + x 7 2 )dp + | ! n^xSdl- + 


|( M s + M d) 


^ 


( 1 ) 


Potential energy V 
° R 

<™\ r 


v = a 


7 R 


0 


(El) T .(u") 2 dr + T(u') 2 dr 

UO 


+ |J o (El) p (v") 2 dp + 


1 

2 


i* 2 

(El) f (x") 2 df + TeQA + K g (v p - v p j 2 


|Kx d o 2 + 


^(^o) 


( 2 ) 


where the dots denote derivatives with respect to time and the primes 
denote derivatives with respect to space or length. In equation (2), 

2 - 

d Q represents rigid-body fuselage translation and p d Q represents rigid- 
body fuselage pitching. 

The chosen deflections are shown in figure 7> along with a, (3, 
and y used as subscripts to designate the following positions on the 
model : 

a intersection of rotor shaft and rotor-disk plane 

3 intersection of rotor shaft and upper pylon frame 

intersection of rotor shaft and horizontal plane of center 
of fuselage 


7 
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With the model mounted in its test configuration, there is a point 
between the front and rear fuselage supports at which a force applied 
vertically produces only fuselage vertical transformation and no fuselage 
longitudinal pitching. The spring constant K x is the force per unit 

deflection for vertical forces applied at this point. The spring con- 
stant Kq is the pitching moment per radian of fuselage rotation about 

this point for pitching moments applied in the longitudinal plane of 
symmetry. The front and rear fuselage supports were mounted so that 
this point coincides with the center of the fuselage (i.e., at f/F = 0.5)* 
The helicopter in free flight would be represented by conditions of K x 

and Kq equal to zero. 


Choice of Modes 

The following equations express the deflections of the model in terms 
of the generalized coordinates and chosen mode shapes: 

For blade first and second flapwise bending, 

u = + a 2 U 2 (3) 

For rotor shaft translation and rotation about the universal joint, 

v = b o v o - <l x 7 (**•) 

For pylon bending and translation, 

W = C]_Wp - px^ (5) 

For fuselage translation, pitching, and bending, 

x = d Q + SqXq + dpxp (6) 

where a-p a 2 , b Q , c-p d D , 3 Q , and d 1 are unknown functions of 
time and u 1 , u 2 , v D , w-j_, x q , and x 1 are the chosen mode shapes. 
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If the inclusion of additional modes in the analysis is desired, 
this is readily accomplished by the inclusion of the appropriate terms 
in equations (3) to (6) for the deflections of the various components. 
For example, third-mode bending of the blades could be included by the 
addition of the term a^u^ l n equation (3)« 

Substitution of equations (3) to (6) into the energy equations 
(eqs. (l) and (2)) results in the following expressions for the kinetic 
and potential energies: 




R 


0 


mr 


a-^u^ + a 2 u 2 + 


1 + 1*0,7 + d i x i, 7 ) 


*1(1 + 1 ^ ~ ^i x i,: 


dr + i 


Q 


mq 


UO 


/. - 2q • 1 

(Vo + 1 Y " d i qx i,: 


/. x - . \ 2 


r p & 

yo + ^o x o ,7 + ^l x l, 7 j 

dqt sJ 

O 


7 2p 


' A2 


A W 1 + d o f- - d lP x l,7 


d o + d o x o, 7 + d l x l, 7 


d P + i 


m f ( d o + d o x o + d l x l) 2<if 


^Sd 


Vo,P + l w l,(3 + 1 y- ■ 2d l Px l ,7 ) 

2l 


V = | 


1 

2 


( d o + l x o,7 + d i x i,r) 

(El) r ^aj_u 1 " + a2U2"^ 2 dr + 


n R 

UO 


Uo 


l?^ a l u l ' + a 2 U 2 dr 


r p pf 

(^ I )p( C l W l") ^dP + \ (El) f ^ d 1 x 1 ") 2 cdf + 


(7) 


^vi + 1 5 ° ' I +K d c i w i,p— o 


b p + 


| K x a o 2 - 4> ^ 


( 8 ) 
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The substitutions 


v a = b o + 5 o ^ " d l C K, 7 
w 3 - V P = c l w l,3 - b o | 


which consist of the fuselage pitching mode, the shaft horizontal 
deflections due to rotation about the universal joint, the total hori- 
zontal deflection of the center of the rotor hub, and the extension of 
the springs and dampers in the upper pylon frame, respectively, are 
used in the derivation of equation (8) . 


Derivation of Frequency Determinant 


Some simplification is obtained by choosing natural modes for the 
uncoupled blade modes up and U 2 and this choice is made in the 

analysis which follows. The differential equations of motion are 
determined by performing the operations indicated by Lagrange's equa- 
tions for free vibrations on the energy equations as follows : 



d/au 

dt Ida 



(9) 


where t denotes time. 


The characteristic equations for free vibrations are found by con- 
sidering harmonic motion, that is, 


ap = ap sin cot 
e.2 = &2 s ^ n ^ * 


do) 


where the amplitudes of ap, & 2 , * * * are denoted by S Q , ap, • • *. 

Substitution of equations (7) and ( 8 ) into equations (9) and making 
use of equations (10) leads to the characteristic equations which are 
given in matrix form as follows: 
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where 


'R 


A-^ = n | dr 

Jo 

P R 

A2 = n m r U2 dr 

Jo 


a 3 = p “pP"! d P 


JO 


A,, = 


0 


ITLfXo 


df 


a 5 = 


Ax = 


Uo 

r*Y 

Jo 


m^Xi df 


mfX 0 X]_ df 


Mi = n 


•R 


p 

ir^Ui dr 


f R 2 

M2 = n m r u 2 dr 

Jo 


m 5 = 


m^x 0 ^df 
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2 df 
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iripP 2 dp 
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The elements of the matrix (eq. (ll)) contain all variables of the 
problem involving distributions of mass and stiffness. Substitution of 
the appropriate values for the model, derived from the model properties 
given in table I and figure 8, yields the following matrix, the elements 
of which contain the natural frequency o>, the rotor speed H, and the 
magnitudes of the generalized coordinates a-^, • • • : 
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Coupled. Frequencies and Mode Shapes 

The characteristic equation of the matrix, which yields the natural 
coupled frequencies, is obtained by setting the determinant of the matrix 
equal to zero. The assignment of unit value to one of the generalized 
coordinates then permits the evaluation of the other coordinates for each 
coupled frequency. The results presented in this report were obtained on 
automatic computing machines by the method of iteration which yields the 
coupled natural frequencies and generalized coordinates simultaneously. 

The mode shapes are then determined by the appropriate combinations of the 
generalized coordinates and the chosen modal functions as follows: 

Fuselage vertical deflections: 

x = + - 2 f) + <*l x l (O) 

Rotor-blade vertical deflections: 

u = S o + ^of 1 - Y 7 ) + a i X l ,7 + *l u l + **2 U 2 ( lk ) 


Rotor longitudinal deflections: 

v a ’ (* 4 ' 3 l Xl ) Q + (15 > 

Pylon longitudinal deflections: 

V = (do| - 2 lXl ^ P + c lWl (16) 

Shaft longitudinal deflections: 

v = (d 0 1 - d x x^ 7 ^q + (17) 

The pylon and shaft vertical deflections are equal to the deflection 
of the fuselage at position 7, that is, x^. 
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For purposes of calculation, the mode shapes used for the components 
Uf, U2> X] _, and w-j_ were derived by fitting suitable approximate mathe- 

matical expressions to the normalized natural uncoupled mode shapes in 
order to obtain analytical expressions which are more readily integrable 
and differentiable. The resulting mode shapes are: 


u n = — - 1.132 sin n— + 0.026 sin 2rt— 
1 R R R 


(18) 


Up = £ - 0.6V7 sin *£ + 0.835 sin 2 n^ + O.I565 sin 3 *- ( 19 ) 

^ R R R R 



(20) 


(21) 


The use of the approximate mode shapes for u^_ and U2 do not 


admit exact conditions of orthogonality for these modes as implied in 
equations (ll) and (12); however, the magnitudes of the cross-coupling 
terms which would exist are negligible and may be ignored. For example 


the coefficient of 



in the elements a^2 and a 21 ^ e ^^rix 


9 


(eq. 12)), if included, would be less than 1 percent of the corresponding 
value for and less than 0.2 percent of the corresponding value for 


£^22 • 


RESULTS AND DISCUSSION 


The following section is divided into two parts. The first part 
presents the results of the calculated frequencies and mode shapes. The 
second part is devoted to a discussion of the comparison between calcu- 
lated and measured results. 


Calculated Results 

Figure 9 shows the uncoupled natural frequencies of the model as a 
function of rotor speed. These frequencies were calculated by assuming 
all elements of the frequency determinant to be zero except those on the 
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principal diagonal and represent physically what would happen if elastic 
deformation occurred only in individual components of the structure. 

For example , the element a^y yields the fuselage uncoupled free-free 

elastic bending frequency under conditions wherein the blades, shaft, 
and pylon are rigid masses. 

The curves of figure 9 show three points where the uncoupled 
natural frequencies coincide. Such a condition, of course, cannot exist 
for a physically coupled system. Figure 10 permits a comparison of the 
calculated uncoupled natural frequencies shown in figure 9 with the 
coupled natural frequencies derived from the complete frequency determi- 
nant. The modal coefficients which correspond to these coupled natural 
frequencies are presented in table II, and the mode shapes for a rotor 
speed of 23 radians per second are shown in the sketch of figure 11. 

The rotor speed of 23 radians per second was chosen for plotting the 
mode shapes because of the substantial coupling indicated by figure 10 
for rotor speeds in this region for the fifth and sixth coupled modes. 

The mode shapes for any other rotor speed are readily obtained by substi- 
tution of the coefficients from table II into equations (13) to (17)* 

The mode shapes shown in figure 11 are normalized by assigning a unit 
value to the point on the model having the maximum deflection for the 
respective mode and computing the relative deflections at other points. 

In plotting the mode shapes, the deflection of only one blade is shown 
since it is assumed in the analysis that all blades behave in a similar 
fashion. The scale for the deflections is substantially exaggerated rela- 
tive to model lengths in order to point out the deflections associated 
with the various modes . The pylon deflections are shown as the dashed 
curves . 

With reference to figures 10 and 11, the principal model deflections 
for the various coupled modes may be summarized as follows: 

Mode 1 - fuselage pitching 

Mode 2 - fuselage translation 

Mode 3 - blade first elastic bending 

Mode b - pylon bending and shaft translation 

Mode 5 - blade second bending and fuselage bending 

Mode 6 - blade second bending, shaft translation, pylon bending, 
and fuselage bending 

Mode 7 - pylon bending 
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Although a seven-degree -of -freedom system was solved, attention is 
focused on the first six modes. The reason for this is as follows: In 

order to obtain a good representation of the lower coupled shaft-pylon 
mode (fourth coupled mode), the flexibility of the pylon was included in 
the analysis, along with the flexibility of the springs in the upper 
pylon frame. As a result, the higher mode, having a coupled frequency of 
approximately 688 radians per second and involving primarily pylon deflec- 
tion, is obtained. This mode is included as a matter of information but 
the mode shape and frequency are only approximate because several 
uncoupled component modes (fuselage second elastic bending, blade third 
and fourth elastic bending, and so forth) having frequencies below the 
uncoupled frequency of the pylon have been neglected since the coupled 
frequencies which would result are above the region of primary concern 
(above the eighth harmonic of the rotor speed) . 


Comparison of Measured and Calculated 
Coupled Natural Frequencies 

The main objective of this investigation was to derive the equations 
of motion for free vibrations by employing logical assumptions as to the 
primary degrees of freedom involved and to compare the calculated natural 
frequencies and mode shapes with those measured on the model to demonstrate 
that sufficient accuracy between calculated and measured results can be 
obtained. 

Both the calculated and measured coupled natural frequencies are 
shown in figure 12 . The calculated frequencies, repeated from figure 10 , 
are shown as the dashed lines, and the measured values are given by the 
symbols. The results presented in figure 12 show good overall agreement 
between calculated and measured values for the coupled natural frequencies 
and indicate that the analytical procedure is adequate and that the modes 
chosen to represent the model structure are the important ones. Thus, it 
is concluded that the coupled natural frequencies and mode shapes can be 
determined by the analytical procedure presented herein with sufficient 
accuracy if the mass and stiffness distributions of the various helicopter 
components are known. 

The analysis did not include damping of the components; however, 
dampers were installed on the model to avoid ground-resonance instabili- 
ties. The results show that substantial damping did exist as evidenced 
by the fact that, for values of the rotor speed where the coupled natural 
frequencies of two modes tended toward each other (e.g.,. modes 5 and 6 at 
2 5 radians per second), it became increasingly difficult to locate the 
individual coupled frequencies on the oscillograph records. Instead, a 
large response of appropriate model components was noted over a rather 
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broad range of frequency of model excitation with but a single peak. 

Such situations are reminiscent of damped coupled systems (page 97 of 
ref . 9) • 

While the coupled natural frequencies were being measured, the 
coupled mode shapes were observed visually by means of stroboscopic 
lighting to ascertain that the relative motions of the components shown 
in figure 11 for the respective modes did exist. For the case of $1 = 0, 
the vibrations of all model components could be readily observed; however, 
when the rotor was turning, rotor-blade motions could not be observed. 
Instead, the blade motions relative to fuselage motions were obtained 
from the oscillograph records of the rotor-blade and fuselage bending 
strains. In all cases, the agreement between the measured and cal- 
culated mode shapes was good but some differences were noted in the loca- 
tion of the node points . 


CONCLUDING REMARKS 


A dynamic model was constructed and tested, and an analytical pro- 
cedure was devised to study the coupled natural frequencies and mode 
shapes of a single-rotor helicopter in the longitudinal plane of sym- 
metry. The analytical study treats a seven-degree-of -freedom system 
involving elastic deformations of the rotor blades, rotor shaft, pylon, 
and fuselage and indicates a straightforward extension for the inclusion 
of additional degrees of freedom if desired. The elements of the fre- 
quency determinant are arranged to permit the uncoupled natural fre- 
quencies and mode shapes, which are often more conveniently determined 
experimentally, to be used to obtain the coupled values. 

A comparison of the uncoupled and coupled natural frequencies for 
the model shows that significant differences exist between these fre- 
quencies for some of the modes, especially at rotor speeds where the 
natural frequencies for two uncoupled modes tend to coincide. 

A comparison of the calculated coupled natural frequencies and mode 
shapes with those measured on the model shows good agreement and suggests 
that the significant variables are properly included in the analysis. 

The coupled natural frequencies and mode shapes can be determined 
by the analytical procedure presented herein with sufficient accuracy if 
the mass and stiffness distributions of the various helicopter components 
are known . 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Field, Va., October 1, 1958* 
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TABLE I.- PARAMETERS OF BASIC CONFIGURATION 


(a) Fuselage parameters 



(b) Blade parameters 





TABLE II.- MODAL COEFFICIENTS 


coefficients for natural coupled frequenc 


1.0000 
5.1120 
• 8334 
.2186 
8.0175 
111.5020 
-19.5985 


1.0000 

.•>776 

- . 0482 
- .0126 
12.1294 
-5. 00 30 
-14. BOl 6 


“(5) 

CD (4) 

<d ( 5) 

“(6) 

“(7) 

ft : 

- 0 




1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

20.3270 

.4114 

20.6682 

-.7166 

-.4435 

2.9330 

4.6829 

-1.56^9 

-.7390 

2.4971 

.7625 

1.1364 

-.3245 

-.1626 

-31.2945 

77.5303 

-.3172 

-.7710 

-.2889 

-.2926 

90.5760 

-.2424 

-.03 ! *3 

.0219 

.2039 

2226.7399 

.1307 

-1.1945 

.4456 

.5967 


ft - 10 radians/sec 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

2.9865 

.3669 

-14.9890 

.3093 

16.6571 

-.8000 

-.4471 

.8302 

-.0427 

-2.4515 

4.7001 

-1.1523 

-.7380 

2.4969 

.2230 

-.0112 

-.6355 

1 . 14 06 

- . 2718 

-.1623 

-31.2919 

8.0264 

12.1871 

-46.8519 

-.314? 

-.6791 

- . 2869 

-.2926 

II3.8165 

-4.6444 | 

-62.3789 

-.2386 

- .6766 

.0255 

.2039 

-38.0329 

-9.9922 

-1393.6470 

.1446 

-.8984 

.4601 

.5976 


ft - 20 radians/sec 


1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1 . 0000 

2.5907 

.4925 

-3.4053 

.1786 

8.8645 

-1.1965 

.8603 

-.0368 

-.6998 

4.7214 

-.9006 

-.7334 

.2264 

-.00963 

-.1799 

1.1459 

-.2080 

-.1608 

8.0482 

12.2524 

-7.7023 

-.3113 

- . 5020 

-.2777 

II7A95I 

-4.2796 

-II.6857 

-.2330 

- . 3640 

.0424 

-22.5308 

-3.0411 

-318.3830 

.1789 

-.2881 

.5177 


ft = 25 radians /sec 

1.0000 1.0000 1.0000 

.4628 -2.6780 .1916 

-.0 396 -0892 4.7258 

-.00955 -.1510 1.1469 

12.2645 -5.5902 -.3104 

-4.2088 -8.6131 -.2315 

-4.1346 -238.4748 .1941 


1.0000 

-.4317 

2.4962 

-31.2838 

-.2924 

.2041 

.6002 


1.0000 

1.0000 

1.0000 

1.0000 

.1316 

6.6173 

-1.4822 

-.4537 

4.7258 

CD 

OJ 

iT\ 

CO 

1 

-.7301 

2.4959 

1.1469 

-.1953 

-.1597 

-31.2802 

- . 3104 

-.4512 

-.2712 

-.2924 

-.2315 

-.2730 

.05429 

.2042 

. 1941 

-.1010 

.5529 

.6014 



TABLE II.- MODAL COEFFICIENTS - Continued 




12.5033 

- 3-9855 

- 1.2631 


- 30.1977 


1.9250 


-.0171 

.6556 


.1*529 

1.4996 



W] 


TABLE II.- MODAL COEFFICIENTS - Concluded 


coefficients for natural coupled frequencies 



1.0000 

1.2713 

.7476 

.2020 

8.1816 

123.232 

- 3.9361 


fi = 30 radians /sec 


1.0000 

.2339 

-.02940 

-.00779 

12.3049 

- 3.9724 

- 1.1273 


1.0000 

1.0000 

-11.9234 

-.4868 

-.6366 

2.4910 

-.1183 

-31.2230 

-.04094 

-.2914 

.4812 

.2037 

1.6266 

.6197 


33 radians/sec 


1.0000 

1.0000 

1.0000 

1.0000 

1 . 0000 

1.0000 

1.0000 

1.1281 

.2064 

.0233 

-.3365 

.3220 

-13.6803 

- .4964 

.7122 

-.0289 

4.3773 

- 1.3437 

-.7307 

-.6231 

2.4896 

.1940 

-.00737 

1.1116 

- . 3680 

-.1660 

-.1090 

-31.2063 

8.2147 

12.3070 

-.3233 

-.8128 

- . 3021 

-.00268 

-.2911 

125.91^ 

- 3.9973 

-.2478 

-.9036 

-.00236 

.3530 

.2061 

-5.0292 

-.9441 

-.6437 

- 27.4793 

.8960 

1.8193 

.6248 


ft - 60 radians /sec 


1.0000 

1.0000 

1.0000 

1 . 0000 

1.0000 

1.0000 

1.0000 

1.0030 

.1827 

.0277 

-.1630 

.2689 

-13.1884 

-.3074 

.6726 

-.0279 

4.6118 

-1.0479 

-.7476 

.6160 

2 . 4880 

. 1891 

-.00734 

1.1200 

-.2436 

- . 1630 

-.0993 

- 31.1879 

8.2310 

12.3084 

-.3177 

-.6043 

-.2933 

.0301 

-.2908 

124.3971 

- 3.9493 

-.2384 

-.9430 

.00983 

.6147 

.2066 

-4 . 3133 

-.8014 

-.3372 

-16.1398 

1 . 2047 

1.9894 

.6303 


70 radians/sec 


1.0000 
.8082 
.9612 
. 1644 
6.3336 
123.9960 


1.0000 
.1431 
- . 0236 
-.00688 
12 . 3097 
- 3.9287 


1.0000 

-17.9897 

-.6034 

-.0779 

.0821 

.7129 


1.0000 

-.3342 

2.4841 

- 31.1431 

-. 290 C 

. 207 C 





Figure 1.- Side view of model showing general methods of shaker installation L-95196.I 

and fuselage mounting. 
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Figure 2.- Front view of the model 


L-95197.1 
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Figure 3 .- Top view of pylon and shaft showing details of spring and damper L-95^98*- 

installations in upper pylon frame • 









Figure 5*- Details of upper pylon and blade lag dampers (assembled and disassembled). L-95202.1 
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Figure 7 • - Sketch showing direction of positive deflections of model test configuration. 
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Top of fuselage 


.116 0463. 


1 Ss , 

f 1 


.114 pounds 

■ ■ i ■ i I - - ■ 1 


0 


5 5.95 7.05 7.6 8.36 in. 


(a) Mass distribution of rotor shaft. 


I / — Top of fuselage 


.414 lb 


r 

\ 

.211 pounds 



0 


5.95 in. 


(b) Mass distribution of pylon. K„ = 84.3 pounds per inch; 
Wp = 584.04 radians per second. 

Figure 8.- Mass distributions of rotor shaft and pylon. 
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Natural 
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Figure 10.- Comparison of calculated coupled and uncoupled natural 
-p _„xes for the basic model configuration. (l-P> euC - 

represent a 1-per ^'^lution frequency, a 2 -per-revolution 
fV - ^ 1 v ■> ^ 



Figure 11.- Calculated mode shapes of coupled system. Rotor speed, 

23 radians per second; pylon deflections are shown as dashed curves 




(d) Mode 4 . l^.kO radians per second. 



(g) Mode 7. 687.9^ radians per second 







0 10 


Figure 12.- Compa] 
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